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This  paper  is  concerned  with  developing  reasonable  criteria  for  the  selection  of  frequence  and 
power  for  very  low  and  low  frequency  (VLET.El  transmitting  stations  The  approach  uses  a  wave 
guide  model  for  low-frequency  propagalion  and  accounts  for  the  variability  of  the  ionosphere  A 
sample  problem  involving  a  hypothetical  transmitter  is  described 


INTRODUCTION 

VLF  LF  radio  signals  are  subject  to  vector  inter¬ 
ference  phenomena  that  can  produce  deep  minima 
in  the  signal  strength  simultaneously  with  rapid 
shifts  of  the  phase  of  the  signal.  These  signal 
minima  are  called  modal  interference  nulls  because 
they  occur  as  the  result  of  destructive  interference 
between  modes  propagating  in  the  Earth-iono- 
sphere  wave  guide.  These  nulls  are  observed  when 
the  signal  amplitude  is  recorded  as  a  function  of 
time,  distance,  or  frequency  within  the  VLF'LF 
band.  These  interference  phenomena  can  be  detri¬ 
mental  to  VLF/LF  communication  systems. 

The  signal  level  varies  with  time  primarily  be¬ 
cause  of  variations  in  the  ionosphere.  Positions  of 
the  modal  interference  minima  also  fluctuate.  The 
position  of  a  null  is  generally  more  stable  during 
normal  daytime  propagation  conditions  than  it  is 
during  the  night.  However,  when  an  ionospheric 
disturbance  occurs,  such  as  those  caused  by  solar 
flares,  the  nulls  may  move  a  100  km  or  may  disap¬ 
pear  completely.  Calculation  of  signal  strength  at  a 
given  location  is  made  difficult  because  it  is  a 
nonlinear  function  of  the  parameters  used  to  de¬ 
scribe  the  Earth-ionosphere  wave  guide.  Hence  it  is 
difficult  to  present  simple  general  conditions  for 
selecting  transmitter  factors  such  as  the  required 
radiated  power.  In  deriving  power  requirements 
which  would  result  in  satisfactory  communications 
with  a  designated  time  availability  at  all  times  of 
day,  it  is  common  to  use  a  worst  case  approach. 
Thus  the  power  requirement  is  based  on  the  condi¬ 
tion  which  results  in  the  lowest  expected  signal-to- 
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noise  ratio  (snr).  Unfortunately,  this  worst  case 
approach  usually  yields  very  high  power  require¬ 
ments.  especially  when  the  effect  of  modal  propa¬ 
gation  nulls  is  considered.  However,  the  deep  nulls 
affect  only  small  geographic  areas  during  the  day¬ 
time.  and  the  impact  on  the  power  required  can  be 
reduced  (derated)  by  neglecting  the  small  percent¬ 
age  of  the  operational  area  affected.  With  this 
technique,  which  will  be  referred  to  as  null  smooth¬ 
ing,  reasonable  radiated  power  requirements  can  be 
derived. 

PROPAGATION  PREDICTION  MODEL 

In  the  propagation  model  employed  in  this  anal¬ 
ysis,  radio  waves  are  considered  to  propagate  in  a 
wave  guide  defined  by  the  space  between  the  Earth 
and  the  ionosphere.  The  model  is  based  on  Budden 
119611  and  was  first  described  by  Puppert  et  al. 
(1967).  In  this  model  both  boundaries  of  the  wave 
guide  may  vary  arbitrarily  over  the  transmission 
path.  These  variations  are  treated  as  a  series  of 
horizontally  homogeneous  segments.  Within  each 
segment  the  ionosphere  is  characterized  by  arbi¬ 
trary  electron  and  ion  density  distributions  and 
collision  frequencies  as  functions  of  height.  The 
lower  boundary  is  considered  to  be  a  smooth  and 
vertically  homogeneous  Earth  which  is  described 
by  its  surface  conductivity  and  dielectric  constant. 
Full  allowance  is  made  for  the  anisotropy  as  a  result 
of  the  Earth’s  geomagnetic  field  and  for  Earth 
curvature.  The  solutions  to  the  mode  equation  are 
obtained  by  numerical  integration  of  the  reflection 
coefficients  through  the  ionosphere.  Field  strength 
is  computed  by  concatenating  the  segments  and 
summing  the  wave  guide  modes  using  a  mode 
conversion  algorithm  \Pappert  and  Snyder .  1972; 
Pappert  and  Morfiit,  1975). 
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The  conductivity  and  dielectric  constant  of  the 
Earth  can  be  specified  with  some  confidence  over 
most  paths  [ International  Radio  Consultative  Com¬ 
mittee  ( CCIR ),  1986].  Likewise,  the  magnitude  and 
direction  of  the  Earth's  magnetic  field  are  well 
known.  The  major  problem  lies  in  the  assignment  of 
parameters  describing  the  ionosphere.  The  distribu¬ 
tion  of  charged  particles  in  the  ionosphere  depends 
in  a  complicated  way  on  latitude,  solar  zenith  angle, 
and  season. 

The  ionospheric  conductivity  parameter  oir  is  a 
function  of  the  charged  particle  density  divided  by 
the  particle-neutral  collision  frequency.  One  of  the 
simplest  ionospheric  profiles  is  an  exponential  vari¬ 
ation  of  conductivity  with  height.  Following  Wait 
and  Spies  [1964],  this  conductivity  is  taken  to  be 

o>r  =  2.5  x  I0'  exp  [filh  —  /i’jJ  (It 

where  h  is  height,  h'  is  the  reference  height,  and  fi 
is  the  gradient  in  inverse  height  units.  The  value  of 
the  electron  density  N  in  electrons  per  cubic  centi¬ 
meter  is  calculated  as  a  function  of  height  as 

V  =  1.43  x  ltf  exp  [pill  -  /i  t  -  0.  l5/i]  (2) 

where  fi  is  in  km  1  and  h'  is  in  kilometers.  The 
collision  frequency  v  in  collisions  per  second  is 
assumed  to  be 

v=  1.816  x  10"  exp  (-0.15/0.  (3) 

The  usefulness  of  this  simple  exponential  iono¬ 
spheric  conductivity  model  for  computing  VLF/LF 
fields  is  demonstrated  by  its  success  in  modeling 
experimentally  measured  data  for  both  daytime  and 
nighttime  propagation  conditions  [Biekel  et  al .. 
1970;  Morfitt  et  al.,  1981].  The  modeling  process 
consists  of  varying  the  values  of  fi  and  h',  calculat¬ 
ing  the  field  strength  as  a  function  of  distance  from 
the  transmitter,  and  comparing  the  calculations 
with  the  measurements.  This  variation  of  iono¬ 
spheric  profile  is  continued  until  acceptable  agree¬ 
ment  between  calculation  and  measurement  is 
reached.  In  general,  the  ionospheric  models  so 
determined  must  be  considered  to  represent  an 
average  ionosphere  since  the  modeling  assumes 
that  the  ionosphere  was  static  during  any  measure¬ 
ment  period.  The  drawback  to  this  approach,  given 
that  the  assumption  of  a  temporarily  static  iono¬ 
sphere  is  valid,  is  that  only  a  single  snapshot  of  the 
ionosphere  results.  Experience  shows  that  this 


1ABLK  I  Optimum  Ionospheric  Parameter  Distribution  for 
Juh 

Parameter  Signal  Statistics 

/  (  I  HO 

At 

tT/j  <M» 

ftV.O 

trl;  U.20 


snapshot  can  change  from  hour  to  hour  and  day  to 
day. 

.SELECT  ION  OF  IONOSPHERIC  PROFILES 

Ferguson  et  ul.  1 1985)  applied  a  statistical  model 
of  the  ionospheric  parameters  fi  and  h'  to  data 
collected  at  Naples.  Italy,  and  La  Maddalena,  Sar¬ 
dinia.  In  this  model,  fi  and  h'  were  assumed  to  be 
distributed  in  a  jointly  normal  distribution.  Best  fit 
values  for  this  model  were  obtained  using  data 
obtained  by  monitoring  the  U.S.  Navy's  LF  station 
located  near  Athens.  The  two  sites  are  on  either 
side  of  a  deep  modal  interference  null  w  hich  occurs 
about  KXX)  km  from  the  transmitter.  The  measured 
signal  strength  mean  and  standard  deviation  were 
used  to  obtain  a  best  fit  mean,  standard  deviation, 
and  correlation  coefficient  for  the  distribution  of  fi 
and  h'.  The  optimum  distribution  of  the  ionospheric 
parameters  yielded  signal  statistics  that  fit  the  mea¬ 
sured  mean  and  standard  deviation  at  both  monitor¬ 
ing  stations  to  within  I  dB.  The  parameters  of  this 
fit  for  July  are  presented  in  Table  I.  where  r  is  the 
correlation  coefficient,  ir ^  is  the  standard  deviation 
of  fi,  and  <iy,  is  the  standard  deviation  of  h' .  The 
parameters  fiu  and  /r(1  are  the  mean  values  of  fi  and 
/t\  respectively. 

Because  the  ionosphere  is  continually  changing 
during  the  day.  the  question  arises  as  to  which 
ionosphere  should  be  used  when  predicting  signal 
strength.  A  common  approach  is  to  use  a  single 
average  ionospheric  profile  to  predict  the  signal 
strength.  The  average  signal  strength  at  a  particular 
location  is  not  always  the  same  as  the  signal  given 
by  the  average  ionosphere.  The  average  ionosphere 
predicts  the  average  signal  strength  quite  well  ex¬ 
cept  in  areas  where  there  are  large  signal  strength 
variations,  such  as  near  modal  interference  nulls. 

To  overcome  some  of  the  weaknesses  associated 
with  using  a  single  average  ionosphere  to  predict 
signal  strength,  multiple  ionospheres  were  used  in 
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this  work  to  include  the  effects  of  variation  of  the 
daytime  ionosphere.  Data  were  generated  for  f3 
from  0.27  to  0.51  km  "  1  in  increments  of  0.02  km  1 
and  for  h’  from  64  to  78  km  in  increments  of  0.5  km. 
Although  the  standard  deviation  of  h‘  in  Table  I 
does  not  require  such  an  extensive  range  of  h\  the 
data  are  useful  for  additional  studies  not  reported 
here.  It  should  be  noted  that  the  volume  of  data 
calculated  for  this  study  is  only  reasonable  because 
the  paths  under  consideration  are  all  seawater  and 
quite  short.  Because  of  the  amount  of  computation 
involved,  it  would  be  very  expensive  to  do  such  a 
study  for  long  paths  over  which  the  ground  conduc¬ 
tivity  changes  frequently. 

PATH  PARAMETERS 

The  application  of  the  propagation  model  to  se¬ 
lection  of  transmitter  frequency  and  power  was 
done  by  generating  calculations  along  two  oppo¬ 
sitely  directed  paths  centered  on  a  hypothetical 
transmitter  site.  The  paths  were  arbitrarily  chosen 
to  form  geographic  bearing  angles  of  96°  and  271°  at 
the  transmitter.  Both  paths  are  taken  to  be  2000  km 
long.  The  snr  is  the  primary  factor  controlling  the 
performance  of  a  communication  system.  The  prop¬ 
agation  prediction  program  used  in  this  study  was 
used  to  calculate  the  signals  along  the  two  paths. 
The  noise  data  were  taken  from  the  atmospheric 
noise  model  described  by  C'CIR  !  1988]  using  July  at 
1600  UT.  The  signal  strength  and  noise  are  assumed 
to  be  separate  Gaussian  distributions  so  that  the  snr 
is  also  Gaussian.  For  purposes  of  discussion  the 
communications  snr  threshold  is  taken  to  be  0  dB, 
Satisfactory  communications  coverage  is  taken  to 
occur  when  this  threshold  is  exceeded  995?  of  the 
time. 

Data  were  generated  from  20  to  60  kHz  in  incre¬ 
ments  of  2  kHz.  The  snr  as  a  function  of  distance  for 
a  radiated  power  of  I  kW  for  these  frequencies  is 
shown  in  Figures  1-4.  The  panels  on  the  left  of  each 
figure  show  results  for  propagation  to  the  east,  and 
those  on  the  right  show  results  for  propagation  to 
the  west.  In  each  panel,  two  curves  are  shown.  The 
data  for  these  curves  were  calculated  using  the 
signal  data  computed  from  the  matrix  of  profiles 
described  above  and  weighted  according  to  the 
Gaussian  distribution  parameters  shown  in  Table  I . 
The  solid  curve  represents  the  mean  snr,  and  the 
dashed  curve  represents  the  snr  which  is  exceeded 
99 %  of  the  time.  There  is  generally  more  than  one 


minimum  in  the  snr  at  each  frequency.  It  is  also 
apparent  that  the  minima  lend  to  move  farther  from 
the  transmitter  as  frequency  increases.  As  the  min¬ 
ima  move  outward,  their  depths  change.  There  is 
also  a  difference  between  eastward  and  westward 
propagation.  It  is  evident  in  Figures  1-4  that  fre¬ 
quency  diversity  can  be  used  to  mitigate  the  effects 
of  nulls  in  critical  areas  by  transmitting  on  frequen¬ 
cies  chosen  such  that  the  positions  of  tne  nulls  from 
the  two  transmitters  complement  each  other. 

CRITERIA  FOR  COVERAGE.  OH  AN  <  IPE.K  A  I  It  IN  AL 
AREA 

ideally,  a  transmitter's  radiated  power  require¬ 
ment  is  based  on  providing  satisfactory  communi¬ 
cations  with  a  designated  time  availability  through¬ 
out  1 005?  of  an  area  of  interest.  However, 
experience  has  shown  that  in  the  regions  of  modal 
interference  nulls  the  signal  level  drops  oft  to  such  a 
low  level  that  very  high  transmitter  power  would  be 
necessary  to  achieve  coverage  of  1005?  of  the  area 
of  interest.  Two  schemes  are  considered:  1 1  >  cover¬ 
age  of  1005?  of  the  area  and  (2)  coverage  using  a 
technique  called  "null  smoothing.”  The  solid  curve 
(999?  time  availability)  in  Figure  5  is  used  to  illus¬ 
trate  how  the  power  requirement  is  determined 
using  these  schemes.  Coverage  of  l(K)r?  of  the  area 
generally  requires  consideration  of  the  snr  at  the 
maximum  distance  and  at  the  locations  of  deep 
minima.  In  Figure  5  the  curve  representing  a  time 
availability  of  995?  has  a  value  of  -15  dB  at  the 
maximum  distance  and  a  value  less  than  -40  dB  in 
the  null  located  near  700  km.  Thus  this  coverage 
criterion  requires  radiation  of  10  MW  in  order  to 
raise  the  snr  in  the  deep  nUI  to  the  communications 
threshold. 

The  second  approach  for  establishing  a  transmit¬ 
ter's  radiated  power  requirement  is  to  use  a  tech¬ 
nique  referred  to  as  "null  smoothing.  "  This  tech¬ 
nique  ignores  the  small  percentage  of  the  area 
affected  by  deep  signal  nulls  when  determining  the 
power  requirement.  Thus  the  minimum  snr  level 
outside  of  the  region  in  which  the  most  severe 
portion  of  the  null  occurs  is  used  as  the  basis  in 
determining  the  power  requirement.  The  compro¬ 
mise  made  using  this  technique  is  that  it  is  accept¬ 
able  for  a  small  percentage  of  the  area  not  to  be 
provided  with  satisfactory  communications  all  of 
the  time.  The  shaded  region  in  Figure  5  corresponds 
to  55f  (100  km)  of  the  total  distance.  As  a  first 
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Fig.  I.  Signal-to-noise  ratio  as  a  function  of  distance  for  ihe  (left -side  panels)  west -east  path  and  the  (right-side 

panels)  east-west  path  for  2IV-28  kff/. 
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Fig.  5.  lllusti.it  ion  of  how  the  power  requirement  iv  determined  with  .md  without  null  smoothing. 


approximation  this  percentage  of  the  path  length  is 
taken  to  represent  a  corresponding  percentage  of 
the  area  which  is  not  covered.  (This  is  not  a  had 
approximation  in  this  case  since  the  nulls  ol  interest 
are  in  the  range  500-1000  km  from  the  transmitter. 
Taking  the  value  to  be  700  km,  as  in  Figure  5.  and 
making  another  crude  approximation  that  the  null 
sweeps  out  a  ring  about  the  transmitter,  the  area  of 
the  null  is  2tt  x  100  x  700.  while  the  corresponding 
total  area  to  be  covered  is  tt  x  2000'.  This  makes 
the  area  of  the  null  about  Wr  of  the  total  area  to  be 
covered).  Ignoring  the  shaded  portion  of  the  snr 
curve  in  the  figure,  the  snr  near  the  null  is  -  13  dB. 
Basing  coverage  on  the  snr  near  the  null  and  at  the 
maximum  distance  results  in  the  power  requirement 
being  determined  by  snr  at  maximum  distance.  The 
power  requirement  becomes  32  kW,  3  orders  of 
magnitude  below  that  required  for  coverage  of 
100%  of  the  area. 

RESULTS  AND  DISCUSSION 

The  power  requirement  as  a  function  of  fre¬ 
quency  for  the  two  coverage  criteria  is  shown  in 


Figure  f>.  This  requirement  is  based  on  maintaining 
the  snr  above  the  threshold  of  0  dB  Wr  of  the  time. 
The  upper  panel  of  Figure  h  is  for  the  path  to  the 
east,  and  the  middle  panel  is  for  the  path  to  the 
west.  In  each  of  these  panels  there  are  five  curves. 
The  dotted  curve  shows  the  power  required  for 
coverage  at  the  maximum  range  (labeled  as  "Max 
Range").  The  data  for  the  other  curves  were  gen¬ 
erated  by  choosing  a  minimum  in  the  snr  at  20  kHz 
and  then  following  i(  from  one  frequency  to  the 
next.  One  pair  of  curves  in  each  panel  was  obtained 
by  following  the  shallow  minimum  near  250  km 
found  in  the  curves  for  20  kHz.  The  small  solid 
circles  connected  by  a  solid  line  (labeled  as  "Min  I 
100%’’)  represent  the  power  required  for  coverage 
of  100%  of  the  area  based  on  this  minimum.  The 
process  of  following  the  minima  was  repeated  with 
application  of  null  smoothing  on  the  same  sequence 
of  minima.  This  third  curve  is  delineated  by  the 
small  solid  circles  enclosed  by  larger  circles  con¬ 
nected  by  a  solid  line  (labeled  as  "Min  I  95 %"). 
Another  pair  of  curves,  show  n  using  crosses  and  a 
dashed  line  (labeled  "Min  2  100%"  and  "Min  2 
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Kig  6.  Power  requirement  as  a  function  of  frequency  for  flop!  Ihe  west  to  east  path,  (middlel  the  east  to  west 

path,  and  (bottom)  combined. 
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95‘ r").  was  generated  starting  with  the  minimum 
found  at  about  600  km  in  the  curves  for  20  kHz. 

The  power  required  for  coverage  of  the  area  at 
each  frequency  is  obtained  by  taking  the  maximum 
of  the  power  required  for  both  directions.  This 
result  is  shown  in  the  bottom  panel  of  Figure  6  in 
w  hich  the  solid  curve  is  for  100#  coverage  and  the 
dashed  curve  is  for  95#  coverage  (based  on  null 
smoothing).  For  I  OCT  coverage  the  pow  er  require¬ 
ment  is  below  100  kW  between  24  and  39  kHz  and 
above  54  kHz.  It  is  below  30  kW  between  26  and  34 
kHz  and  above  56  kHz.  The  minimum  power  re¬ 
quirement  is  22  kW  at  32  kHz.  Coverage  of  95#  of 
the  area  is  obtained  with  a  power  requirement 
below  50  kW  for  all  frequencies  above  25  kHz  with 
a  minimum  of  Ik  kW  at  38  kHz. 

The  selection  of  frequency  and  power  for  a  trans¬ 
mitting  station  is  a  complicated  task  involving 
trade-otfs  between  the  engineering  and  cost  factors 
and  the  requirements  for  acceptable  communica¬ 
tions  and  lime  availability.  Designing  to  a  worst 
case  is  usually  unsatisfactory  in  the  light  of  the 
uncertainties  in  signal  variability  at  low  frequen¬ 
cies.  A  more  reasonable  approach  is  to  accept 
slightly  less  than  HHY'r  coverage  of  an  area.  It 
remains  to  evaluate  by  process  of  measurement  and 
modeling  the  overall  impact  of  this  approach. 
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